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'  Abstract 

The  problem  of  the  inversion  of  irradiance  measurements  for  the  inherent  optical  properties  of 
a  hydrosol  such  as  the  world  ocean  has  been  examined  with  the  NOARL  optical  model.  This 
model  is  a  Monte  Carlo  simulation  of  the  radiative  transfer  equation  that,  given  the  optical 
properties  of  a  hydrosol,  generates  the  zonal  radiances  and  the  commonly  measured  irradiances 
of  that  hydrosol.  This  information  allows  us  to  test  the  optical  inversion  capability  of  any  irradiance 
model.  We  have  demonstrated  that  the  two-flow  model,  used  often  in  atmospheric  optics,  cannot 
be  inverted  to  yield  the  inherent  backscatter  coefficient  of  natural  hydrosols  due  to  the  extreme 
asymmetry  of  the  hydrosol  volume  scattering  function.  The  asymmetry  of  the  volume  scattering 
function  interacting  with  the  radiance  distribution  introduces  parameters  called  shape  factors  into 
the  models  for  irradiance  inversion.  The  presence  of  shape  factors  in  an  irradiance  model  renders 
it  unsolvable;  thus,  simplifying  assumptions  are  required  for  an  inverse  solution-We  inquire -tnw 
the  possibilities  for  irradiance  inversion  given  the  existence  of  shape  factors. 

/ 
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There  is  great  interest  in  the  problem  of 
inversion  of  irradiance  measurements  in  a 
hydrosol  (pure  water  plus  solutes  plus  sus- 
pensoids).  Is  it  possible,  given  a  suite  of 
irradiance  measurements  of  the  submarine 
light  field,  to  invert  these  measurements  to 
determine  the  optical  properties  of  a  hy¬ 
drosol  such  as  the  world  ocean?  The  optical 
properties  of  in  terest  are  the  absorption  coef¬ 
ficient  and  the  backscatter  coefficient  of  the 
hydrosol.  In  many  cases,  the  primary  con¬ 
tribution  to  absorption  and  scattering  comes 
from  the  suspended  particles,  which  are  pri¬ 
marily  living  cells.  It  is  often  possible  to 
differentiate  the  living  from  the  nonliving 
components  by  the  variation  in  these  op¬ 
tical  properties  at  different  wavelengths.  The 
absorption  and  backscatter  coefficients  are 
linear  functions  of  the  hydrosol  components 
they  represent.  Thus,  they  offer  great  poten- 
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rial  for  theoretical  prediction  of  potential 
productivity  of  different  oceanic  regions, 
ecological  energy  budgets,  and  meteorolog¬ 
ical  energy  budgets. 

The  problem  of  optical  inversion  begins 
with  consideration  of  the  radiative  transfer 
equation,  accepted  as  the  standard  mathe¬ 
matical  description  of  atmospheric  and  sub¬ 
marine  light  fields  (Chandrasekhar  1960; 
Preisendorfer  1976).  However,  this  equa¬ 
tion  has  no  simple  analytical  solution  that 
is  applicable  to  hydrosols  and  requires  a  large 
number  of  light  radiance  measurements  to 
solve  it  for  the  inversion  problem  (Zaneveld 
1 974).  These  measurements  are  difficult  and 
expensive  to  obtain,  although  recent  de¬ 
velopments  in  microprocessor-controlled 
instrumentation  may  allow  the  measure¬ 
ments  required  for  the  radiative  transfer 
equation  in  the  future  (Voss  1988).  If,  how¬ 
ever,  the  optical  coefficients  of  a  hydrosol 
are  known,  the  radiant  fluxes  predicted  from 
the  radiative  transfer  equation  can  be  sim¬ 
ulated  with  Monte  Carlo  models  such  as  the 
NOARL  optical  model.  Thus  any  algo¬ 
rithm  that  inverts  radiant  flux  to  estimate 
optical  parameters  can  be  independently 
tested  with  this  Monte  Carlo  model. 

An  approach  to  the  problem  of  inverting 
radiant  fluxes  to  obtain  optical  coefficients 
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has  been  made  by  integrating  the  radiative 
transfer  equation.  This  yields  simpler  equa¬ 
tions  in  terms  of  the  irradiances  that  can  be 
conveniently  measured  in  the  field.  The  in¬ 
tegration  of  the  radiative  transfer  equation 
over  all  viewing  angles  yields  the  Gershun 
(1939)  equation.  Integrating  the  radiative 
transfer  equation  over  the  downwelling 
hemisphere  and  over  the  upwelling  hemi¬ 
sphere  of  three-dimensional  space  yields 
Schuster-type  equations,  originally  pro¬ 
posed  by  Arthur  Schuster  for  stellar  atmo¬ 
spheres  in  1905.  The  Schuster  equations 
were  the  original  stimulus  to  the  develop¬ 
ment  of  the  radiative  transfer  equation  and 
are  still  used  to  study  the  problem  of  optical 
inversion  in  aerosols  and  hydrosols.  The 
most  recently  proposed  integrations  of  these 
equations  for  hydrosols  are  from  Aas  (1987) 
and  Preisendorfer  and  Mobley  (1984).  Za- 
neveld  (1982)  has  proposed  integrations  of 
the  radiative  transfer  equation  that  yield  an 
equation  utilizing  the  irradiances  of  the  light 
field  and  the  upwelling  radiance. 

In  this  paper  we  investigate  the  possibil¬ 
ities  for  inversion  of  simple  irradiance  mea¬ 
surements  in  the  field  to  obtain  the  optical 
properties  of  Case  1  waters,  which  comprise 
98%  of  the  world  ocean  (Morel  1988). 

Optical  properties  and  the 
radiative  transfer  equation 

All  optical  properties  and  radiant  flux 
quantities  discussed  here  are  assumed 
wavelength-specific  so  that  wavelength  is  not 
delineated  in  the  properties  and  their  deri¬ 
vations.  Definitions  of  the  radiance  fluxes 
measured  in  the  marine  hydrosol  are  given 
in  the  list  of  symbols. 

The  optical  properties  we  wish  to  obtain 
from  the  various  inversion  models  depend 
only  on  the  molecular  composition  and 
structure  of  the  constituents  of  a  hydrosol 
and  are  unaffected  by  changes  in  the  distri¬ 
bution  of  radiant  flux.  Such  optical  prop¬ 
erties  were  termed  inherent  optical  prop¬ 
erties  by  Preisendorfer  (1976)  because  they 
are  inherent  to  the  medium  being  studied. 
These  optical  coefficients  are  strictly  linear¬ 
ly  related  to  the  concentration  of  the  various 
components  of  the  hydrosol  and  their  effects 
on  light  field  transmission  are  thus  precisely 
additive.  The  inherent  optical  properties  on 


which  we  will  focus  are  the  absorption  coef¬ 
ficient  a  and  the  backscatter  coefficient  bb 
while  the  derivations  will  require  consid¬ 
eration  of  the  volume  scattering  function 

m. 

We  will  also  consider  the  optical  prop¬ 
erties  of  the  radiation  field  of  the  hydrosol. 
Among  the  most  important  of  these  prop¬ 
erties  are  the  mean  cosines  (Schellenberger 
1965)  of  variously  defined  average  photon 
paths.  The  mean  cosines  are  used  in  the 
inversion  models  and  are  functions  of  the 
distribution  of  radiant  flux  in  the  marine 
hydrosol.  These  properties  are  descriptive 
of  the  structure  of  the  light  field  penetrating 
the  hydrosol  and  have  thus  been  termed  the 
light  field  properties  by  Kirk  (1983). 

Finally,  there  is  a  class  of  optical  prop¬ 
erties  that  represents  combinations  of  in¬ 
herent  optical  properties  and  light  field 
structure.  Preisendorfer  (1976)  termed  this 
combination  apparent  optical  properties.  We 
will  investigate  the  shape  factors  which,  be¬ 
cause  they  are  a  function  of  the  hydrosol 
volume  scattering  function  and  the  radiance 
distribution,  are  apparent  optical  proper¬ 
ties.  These  parameters  are  derived  below. 
Also  in  this  class  of  optical  properties  are 
the  radiance  backscatter  coefficients  given 
in  the  list  of  symbols. 

Radiative  transfer  equation— This  equa¬ 
tion  is  an  integro-diflerential  description  of 
the  absorption  and  scattering  of  the  radiant 
energy  in  a  specified  radiance  at  a  point  z 
in  the  hydrosol 


dL(z;  6,  <f>) 
dz 


cos  6  =  -cL(z\  6 ,  <t>) 


+  L*{z ;  6,  4>) 


(1) 


where 
L\z-  6 ,  <)>) 


0(z;  6,  0,  0\  4>')L(z ;  O',  do.'. 


z  is  the  depth  in  meters,  6  the  zenith  angle, 
<f>  the  azimuth  angle,  c  =  a  +  h,  and  d«  * 
sin  6  dd  d <t>.  The  first  term  on  the  right-hand 
side  of  Eq.  1  is  the  loss  of  radiant  energy 
from  the  radiance  L(z;  0,  <t>)  due  to  absorp¬ 
tion  and  scattering.  The  second  term— the 
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L(r,  8,  <j>) 


EJLz),  E.(z) 


E(z) 


E  0 d(z  ). 

E„  Jz) 


Eo(z) 


»/z),  Uz) 


Significant  symbols 

Radiance— radiant  fiux  per  unit  area  per 
unit  solid  angle  received  at  depth  z  from 
the  bearing  of  zenith  angle  8  and  azimuth 
angle  p  by  a  tube  or  lens  system  pointing 
toward  these  coordinates,  W  m-2  sr"' 


Down  welling  and  upwellingirradiance—  ..  . 
downward  or  upward  flow  or  stream  of 
radiant  energy  at  depth  z  which  repre¬ 
sents  the  summed  vertical  downward-  or 
upward-moving  components  of  the  ra¬ 
diances 


L(z\  8,  p) cos  8  du,  W  m"2 


J'  L(z:  8,  p) cos  8  du,  W  m'2 
2» 


a 


Downwelling  vector  irradiance— the  net  j 
flow  of  radiant  energy  at  depth  z\  the  up- 
welling  irradiance  subracled  from  the 
downwelling  irradiance 


Hz',  8,  ip) cos  8  du  =  E4  -  W  m  2 


4« 


K 


Downwelling  and  upwelling  scalar  irra- 
diance— the  summation  of  all  radiances 
in  the  downwelling  or  upwelling  hemi¬ 
spheres  at  depth  r 


■; 


L(z\  8.  ip)  du,  W  m  2 


b, 


J  L(r;  #,  p)  du,  W  m  2 


m 


Scalar  irradiance  — the  summation  of  all 
radiances  of  three-dimensional  space  at 
depth  z 


L(z\  8 ,  ip)  du,  W  m  2 


Mean  cosine  of  the  average  path  traveled 
by  downward-  or  upward-moving  pho¬ 
tons  at  depth  z 


f  E(z;  8, 


<p)cos  8  du 


f  L 


8 ,  p)  du 


E/z) 

EaAz) 


I 


L(z ;  8,  <p) cos  8  du 


2wy 

J~  L(z\  8,  <p)  du 


EJjz) 

Eou(z) 


Mean  cosine  of  the  average  path  traveled 
by  all  photons  at  depth  z 


I  L(z\  8,  p)c os  8  du 

i _ _  K(z) 

J  L(z\  8,  p)  du  £°(Z) 


Absorption  coefficient— a  measure  of  ra¬ 
diant  energy  absorbed  from  a  beam  of 
radiant  energy  and  converted  into  another 
form  by  an  infinitesimally  thin  layer  of 
hydrosol,  nr1 

Total  scattering  coefficient — a  measure  of 
radiant  energy  diverted  from  a  beam  of 
radiant  energy  by  an  infinitesimally  thin 
layer  of  hydrosol 

=  bh  +  bf,  m"1 

Backseat  ter  coefficient— a  measure  of  ra¬ 
diant  energy  diverted  by  an  infinitesi¬ 
mally  thin  layer  of  hydrosol  from  a  ref¬ 
erence  beam  of  radiant  energy  and  moving 
in  a  direction  90°-180°  to  the  reference 
beam  direction,  m  1 

Forward  scattering  coefficient— a  mea¬ 
sure  of  radiant  energy  diverted  from  a 
reference  beam  by  an  infinitesimal  thick¬ 
ness  of  hydrosol  and  moving  in  a  direc¬ 
tion  0°-90°  to  the  reference  beam  direc¬ 
tion,  m_l 

Volume  scattering  function— a  measure 
of  radiant  intensity  diverted  from  a  beam 
of  radiant  energy  by  an  elemental  volume 
of  hydrosol  into  the  direction  indicated 
by  the  angle  t  relative  to  the  orientation 
of  the  incident  beam,  m  1  sr1 

Downwelling  and  upwelling  radiance 
backscatter  coefficients— measures  of  the 
amount  o T  radiant  energy  backscattered 
from  the  downwelling  or  upwelling  irra¬ 
diance  by  an  infinitesimal  thickness  of 
the  medium,  rrr1 

Downwelling  and  upwelling  shape  fac¬ 
tors— the  normalized  contributions  to  the 
downwelling  or  upwelling  radiance  back¬ 
scatter  coefficient  from  the  forward  scat¬ 
tering  lobe  of  the  hydrosol  volume  scat- 
tering  function _ 
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path  function— is  the  gain  of  radiant  energy 
into  the  radiance  L(z;  8,  <t>)  due  to  scattering 
of  radiant  energy  out  of  all  the  other  radi¬ 
ances  in  three-dimensional  space.  The  ze¬ 
nith  angle  8'  and  the  azimuth  angle  <j>'  refer 
to  the  orientations  in  space  of  the  contrib¬ 
uting  radiances.  The  notation  for  the  vol¬ 
ume  scattering  function  differs  from  the  list 
of  symbols  in  that  the  angle  ^  of  intensity 
scattered  into  the  radiance  L(z\  8,  <f>)  from 
radiance  L(z;  S',  </>')  is  the  angle  between  the 
spatial  coordinates  8,  <f>  and  6',  <f>'. 

Integrated  forms  of  the  radiative  transfer 
equation  for  irradiances—  Because  of  the  dif¬ 
ficulties  in  inverting  the  full  radiative  trans¬ 
fer  equation,  various  integrations  of  this 
equation  have  been  proposed  that  use  the 
irradiances  defined  in  the  list  of  symbols. 
Integrating  Eq.  1  over  all  viewing  angles  of 
three-dimensional  space,  we  obtain 

dE(z) 

~~T  =  -a(z)E0(z),  (2) 

which  is  Gershun’s  (1939)  equation— an 
expression  of  the  net  flow  of  radiant  energy 
in  a  hydrosol  and  a  fundamental  physical 
statement  of  conservation  of  energy  in  the 
divergence  of  a  radiative  field.  The  inver¬ 
sion  of  this  equation  with  the  three-param¬ 
eter  model  (Stavn  1981,  1987)  gives  only 
the  absorption  coefficient.  Additional  inte¬ 
grations  are  necessary  to  get  expressions  in 
terms  of  the  scattering  coefficients. 

Integrating  Eq.  1  first  over  the  down- 
welling  hemisphere  and  then  separately  over 
the  upwelling  hemisphere  yields  two  differ¬ 
ential  equations  in  terms  of  the  downwelling 
and  the  upwelling  streams  of  radiant  flux 


d EJ& 

d  z 


a(z) 

£,y(z) 


EJ& 


+ 


Hz) 


EJ,z) 


b(z) 

Uz) 


EAz) 


(3) 


d Eu(z) 
d  z 


a(z) 

Mz) 


Eu(z)  - 


Hz) 

M„(z) 


Eu(z) 


+  ^^z) 

v-Az) 


(4) 


where  B(z)  is  the  mean  backscatter  coeffi¬ 
cient  at  depth  z.  Equations  3  and  4  are  the 


integrations  proposed  by  Preisendorfer  and 
Mobley  (1984),  which  can  be  inverted  and 
solved  as  a  system  of  two  equations  in  two 
unknowns.  They  represent  the  first  real  ad¬ 
vance  in  the  application  of  the  concept  of 
the  two-flow  equations  of  Schuster  (1905) 
to  hydrosols.  Previous  attempts  combined 
the  light  field  properties,  inherent  optical 
properties,  and  apparent  optical  properties 
in  the  right-hand  side  of  Eq.  3  and  4  to  yield 
poorly  defined  optical  coefficients  that  could 
not  be  applied  in  any  general  sense  to  any 
part  of  the  ocean. 

The  Preisendorfer  and  Mobley  equations 
are  stated  in  terms  of  a  true  hydrosol  ab¬ 
sorption  coefficient  and  a  mean  backscatter 
coefficient  applied  to  both  irradiance 
streams.  They  defined  the  mean  backscatter 
coefficient  as  a  close  approximation  of  the 
hydrosol  backscatter  coefficient.  The  back¬ 
scatter  coefficient  of  the  hydrosol  has  been 
shown  by  Gordon  et  al.  (1975)  to  be  the 
most  important  aspect  of  the  hydrosol  scat¬ 
tering  coefficient  in  terms  of  altering  the  ra¬ 
diance  distribution  and  therefore  the  nature 
of  the  radiant  flux  in  the  marine  hydrosol. 

The  mean  backscatter  coefficient  approx¬ 
imates  the  hydrosol  backscatter  coefficient 
when  the  volume  scattering  function  of  the 
hydrosol  is  uniform  or  symmetrical.  The 
essential  justification  for  the  mean  back¬ 
scatter  coefficient  approximation  is  the 
nearly  uniform  shape  of  the  backscattering 
lobe  of  the  hydrosol  volume  scattering  func¬ 
tion.  This  allows  the  backscattering  of  ra¬ 
diant  flux  to  be  treated  in  the  same  way  as 
uniform  scattering,  i.e.  we  are  dealing  with 
half  of  a  uniform  scattering  function.  How¬ 
ever,  a  problem  arises  because  of  the  large 
forward  scattering  lobe  of  particles  sus¬ 
pended  in  natural  waters  (Petzold  1 972;  Jer- 
lov  1976)  and  the  shape  of  the  radiance  dis¬ 
tribution.  It  is  well  known  (Smith  1974; 
Jerlov  1976)  that  the  radiance  distribution 
of  the  downwelling  hemisphere  is  elongate 
and  the  maximum  radiances  occur  within 
40°  of  the  vertical  axis.  In  contrast,  the  max¬ 
imum  of  the  radiance  distribution  in  the 
upwelling  hemisphere  occurs  near  the  hor¬ 
izontal  plane.  The  backscattering  of  the 
downwelling  and  upwelling  irradiances  aris¬ 
es  not  only  from  the  backscattering  lobe  of 
the  inherent  volume  scattering  function  of 
the  hydrosol  but  also  from  the  forward  scat- 
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tering  lobe.  The  mean  backscatter  coeffi¬ 
cient,  strictly  speaking,  is  therefore  not  di¬ 
rectly  referable  to  the  backscattering  lobe  of 
the  inherent  volume  scattering  function  of 
the  hydrosol  (Fig.  1).  Preisendorfer  and 
Mobley  were  well  aware  of  this  and  sug¬ 
gested  that  how  far  a  single  mean  backscat¬ 
ter  coefficient  for  the  two  irradiance  streams 
deviates  from  being  an  inherent  optical 
property,  an  approximation  of  the  hydrosol 
backscatter  coefficient,  and  a  quantitative 
estimate  of  this  coefficient  was  worthy  of 
theoretical  and  experimental  investigation. 

Aas  (1987)  also  integrated  the  radiative 
transfer  equation  over  the  downwelling  and 
upwelling  hemispheres  to  yield  the  follow¬ 
ing  pair  of  equations 


be  inverted  because  it  is  a  statement  of  two 
equations  in  four  unknowns. 

Defining  the  radiance  backscatter  coeffi¬ 
cients  in  terms  of  Preisendorfer  and  Mob¬ 
ley’s  notation,  we  review  the  particular  in¬ 
tegrations  of  the  radiative  transfer  equation 
that  produce  them 


bd  =  r Job 
i>u  =  rj)h 


and 

bAz) 


■(III. 


0(0,  <t>,  0',  <t>’)  do; 


•L(z;  0',  <(>')  da>'f  /  EoAz) 


(7) 

(8) 


(9) 


AEAz) 

dr 


a(z) 

vA?) 


EA?) 


rAz)bh{z) 

aAz) 


EAz) 


+ 


^u(^)^fc(~)  r’ 

- J  ...  Cut-) 

Mu(-) 


(5) 


and 


d/:„(r) 

dr 


a(  r) 


EAs) 


+ 


Mu(2) 
r„{:)bh(:) 
Au(^) 

rj,z)bh(z) 

fiAz) 


Eu(z) 

EA=) 


(6) 


where  bh  is  the  hydrosol  backscatter  coeffi¬ 
cient  and  rd  and  r„  are  the  shape  factors. 
The  combined  coefficients  (rjy,,)  nd  (rjbh) 
can  be  called  the  radiance  backsc.ter  coef¬ 
ficients  which  indicate  the  amount  of  irra¬ 
diance  backscattered  from  each  irradiance 
stream  divided  by  an  infinitesimal  thickness 
of  the  hydrosol.  These  coefficients  serve  the 
same  function  in  Aas’  integration  as  the 
mean  backscatter  coefficient  in  the  Preisen¬ 
dorfer  and  Mobley  integration  and  thus  em¬ 
phasize  the  fact  that  certain  assumptions 
have  to  be  made  to  obtain  a  mean  back- 
scatter  coefficient.  Aas’  integration  cannot 


{-(z)  ■  II. 

•L(r;0\*')d«'|/  E0Az)  (10) 

where  2vd  and  2*-,,  indicate  integrations  of 
the  three-dimensional  space  over  the  down- 
welling  and  upwelling  hemispheres.  From 
the  above  integrations  we  see  that  the  back¬ 
scatter  coefficients  of  the  irradiance  streams 
involve  truncated  portions  of  the  entire  hy¬ 
drosol  volume  scattering  function.  There¬ 
fore  the  forward  lobe  of  the  volume  scat¬ 
tering  function  enters  into  the  radiance 
backscatter  coefficients.  The  shape  factor  is 
then  seen  as  a  measure  of  the  contribution 
of  the  forward  scattering  lobe  of  the  hydro- 
sol,  for  a  given  radiance  distribution,  to  the 
radiance  backscatter  coefficient.  The  shape 
factor  normalizes  the  contribution  of  the 
forward  scattering  lobe  of  the  hydrosol  to 
the  backscattering  lobe.  This  results  in  a 
mathematical  expression  that  ultimately  re¬ 
lates  changes  in  the  irradiances  to  the  back¬ 
scatter  coefficient  of  the  hydrosol. 

The  inversion  proposed  by  Preisendorfer 
and  Mobley  for  obtaining  the  backscatter 
coefficient  of  the  hydrosol  is  possible  only 
if  the  mean  backscatter  coefficient  defined 
from  the  two  irradiance  fluxes  of  the  two- 
flow  model  has  the  following  properties  rel¬ 
ative  to  the  integration  proposed  by  Aas: 


[J  . 
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Contributions  to  Upweiling  Backscattered 
Irradiance  from  Radiance  LfoB'^') 


From  Forward  From  Backscatterlng  Lobe 

Scattering  Lobe 


Fig.  1.  Downward-moving  radiance  at  depth  2,  L(z\  6',  <p'),  undergoing  scattering  by  the  hydrosol.  Volume 
scattering  function  of  intensities  diverted  from  the  radiance  in  plane  of  paper  indicated  by  envelope  of  the  thm 
line.  Scattered  intensities  contributing  to  backscattered  radiances  come  from  shaded  portion  of  volume  scattering 
function.  Contribution  of  £(r;  8',  <t>)  to  two  upward-moving  radiances,  L(z\  8„  0,)  and  L(z:  $2,  02),  indicated  by 
white  lines  in  the  shaded  portion  of  volume  scattering  function.  Backscattered  radiances  contributing  to  upweiling 
irradiance  at  depth  2  dominated  by  forward  lobe  of  volume  scattering  function. 


h  ~  Bh  =  hu, 

which  is  true  only  when 
rd  =  ru. 

However,  for  the  mean  backscatter  coef¬ 
ficient  to  be  equal  to  the  inherent  backscat¬ 
ter  coefficient  of  the  hydrosol,  the  following 
must  be  true 

rd  =  ru  =  1 .0. 

That  is,  the  mean  backscatter  coefficient  is 
a  valid  estimate  of  the  hydrosol  backscatter 
coefficient  only  when  the  foward  scattering 
lobe  of  the  hydrosol  does  not  make  a  sig¬ 
nificant  contribution  to  the  backscattered 
flux  of  the  irradiance  field.  Aas  thus  pro¬ 
vided  the  analysis  that  allows  the  question 
of  the  inversion  of  two-flow-type  models  to 
be  posed  very  simply  and  quantitatively: 
how  much  do  the  shape  factors  (rd  and  rj 


deviate  from  1 .0  in  ocean  waters?  His  initial 
calculations  of  rd  and  r„  indicated  that  in 
the  most  turbid  coastal  waters  (Case  2 
waters)  the  Preisendorfer  and  Mobley  so¬ 
lution  was  probably  not  valid.  Aas  specu¬ 
lated  that  in  the  clearest  ocean  waters  the 
deviation  from  1 .0  might  not  be  too  great 
and  the  two-flow  model  could  be  inverted. 
Based  on  the  work  of  Aas,  we  extend  the 
analysis  of  rd  and  ru  to  Case  1  waters  with 
the  NOARL  optical  model.  This  simula¬ 
tion  of  the  radiative  transfer  equation  al¬ 
lows  a  quantitative  calculation  of  the  mag¬ 
nitude  of  the  shape  factor  and  of  the  relative 
contributions  of  the  various  components  of 
the  marine  hydrosol,  a  prediction  of  effects 
based  on  known  optical  properties  of  the 
open  ocean  as  reported  extensively  by  Mo¬ 
rel  and  coworkers,  and  inclusion  of  the  effect 
of  hydrosol  absorption  which  always  re¬ 
duces  or  ameliorates  the  effects  of  scatter- 
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ing.  With  this  approach  we  can  investigate 
the  question  of  optical  inversion  in  natural 
hydrosols  and  determine  the  limitations  of 
the  simple  irradiance  solutions  so  far  pro¬ 
posed  for  the  radiative  transfer  equation. 

Case  1  waters:  Optical  parameters 
for  the  NOARL  Monte  Carlo  simulation 

Morel’s  optical  classification  of  ocean 
waters  represents  a  synthesis  of  various  as¬ 
pects  of  the  remarkably  useful  and  fruitful 
hypothesis  on  the  biological  origin  of  the 
optical  properties  of  Case  1  ocean  waters. 
The  primary  sources  of  this  hypothesis  are 
Smith  and  Baker  (1978),  Morel  and  Prieur 
(1977),  and  Morel  (1980).  The  aspect  of  the 
Morel  optical  classification  of  ocean  waters 
that  we  use  here  is  the  association  of  the 
total  scattering  coefficient  of  the  hydrosol  at 
550  nm  with  the  concentration  of  the  chlo¬ 
rophyll-like  pigments  of  the  hydrosol,  after 
the  suggestion  of  Smith  and  Baker  about  the 
importance  of  chlorophyll-like  pigments  in 
Case  1  waters.  The  wavelength  chosen  for 
the  simulations  was  440  nm,  which  is  close 
to  the  peak  absorption  for  Chi  a.  This  wave¬ 
length  has  been  used  extensively  for  ocean 
optical  modeling. 

Scattering  parameters  of  the  simula¬ 
tions—  All  of  our  simulations  begin  with  the 
restricted  version  of  the  NOARL  blue  water 
model  (Stavn  and  Weidemann  1 988a)  with¬ 
out  emission  sources.  This  model  is  equiv¬ 
alent  to  the  Morel  blue  water  model  type  3 
(Morel  and  Prieur  1977).  The  basis  of  the 
model  is,  of  course,  the  optical  properties 
of  molecular  water;  the  total  scattering  coef¬ 
ficient  and  volume  scattering  function  of  the 
water  molecule  (determined  from  fluctua¬ 
tion  theory)  are  taken  from  Morel  (1974). 
The  other  component  of  the  model  is  finely 
divided  quartzlike  material  which  presum¬ 
ably  has  an  aeolian  source  (Brown  and  Gor¬ 
don  1973).  The  total  scattering  coefficient 
and  volume  scattering  coefficient  for  this 
material  are  from  Kullenberg  (1968)  and 
Gordon  et  al.  (1975). 

For  the  living  component  of  the  simula¬ 
tion  we  chose  the  alga  Platy monos  sp.;  its 
optical  properties  are  taken  from  Bricaud 
and  Morel  ( 1 986).  The  total  scattering  coef¬ 
ficient  and  volume  scattering  function  were 


then  determined  from  these  properties.  The 
total  scattering  coefficient  was  normalized 
to  the  concentration  of  chlorophyll-like  pig¬ 
ments.  The  cell  size  distribution  from  Bri¬ 
caud  and  Morel  had  a  median  value  of  4 
jtm  in  diameter.  Multispecific  assemblages 
can  be  modeled  for  the  generation  of  Case 
1  waters  as  the  optical  properties  of  other 
species  become  available. 

The  total  scattering  coefficient  for  organic 
detritus  was  estimated  from  the  optical 
model  of  Prieur  and  Sathyendranath  (1981). 
Their  model  scattering  coefficients  are  nor¬ 
malized  to  550  nm  and  converted  to  other 
wavelengths  from  the  approximately  X-1  re¬ 
lationship  of  Morel  (Morel  1973;  Morel  and 
Prieur  1977).  Prieur  and  Sathyendranath  es¬ 
timated  the  quantitative  contribution  of  or¬ 
ganic  detritus  to  the  hydrosol  total  scatter¬ 
ing  coefficient  by  assuming  the  following. 
The  minimal  observed  hydrosol  scattering 
coefficients  at  550  nm  for  all  Case  1  waters 
are  due  to  living  algal  cells  plus  water  and 
quartz.  The  maximum  hydrosol  scattering 
coefficients  observed  for  these  waters  result 
from  the  contribution  of  suspended  organic 
detritus  and  bacteria  added  to  the  algal  cell 
component.  Such  estimates  are  based  on  the 
empirical  limits  established  by  Morel  (Mo¬ 
rel  1980;  Gordon  and  Morel  1983)  and  il¬ 
lustrated  in  Fig.  2.  The  difference  between 
these  two  values  constitutes  the  total  scat¬ 
tering  coefficient  for  suspended  organic  de¬ 
tritus. 

We  have  modified  this  procedure  as  fol¬ 
lows.  For  a  given  pigment  concentration  the 
total  scattering  coefficient  at  550  nm  due  to 
Platymonas  was  calculated  from  Bricaud 
and  Morel  ( 1 986)  and  then  subtracted  from 
the  maximum  hydrosol  total  scattering  coef¬ 
ficient  at  550  nm  with  water  and  quartz  scat¬ 
ter  removed.  The  result  is  an  estimate  of 
the  total  scattering  coefficient  of  particulate 
organic  detritus.  This  coefficient  was  then 
converted  to  its  ,ralue  at  440  nm.  The  hy¬ 
drosol  total  scattering  coefficient  calculated 
for  Platymonas  plus  the  blue  water  com¬ 
ponent  exceeded  Morel’s  upper  empirical 
limit  at  higher  chlorophyll  concentrations. 
Accordingly,  an  empirical  curve  fit  between 
the  organic  detritus  total  scattering  coeffi¬ 
cient  and  the  chlorophyll-like  pigments  at 
lower  concentrations  allowed  extrapolation 
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Chi  +  Pheo  mg/m3 

Fig.  2.  Total  scattering  coefficients  for  hydrosots  of  NOARL  simulations.  Scattering  coefficients  for  550  nm 
shown  to  facilitate  comparison  with  the  Morel  optical  classification  of  ocean  waters.  Dashed  line  ML  represents 
Morel’s  lower  limit  for  total  scattering  coefficient  predicted  from  chlorophyll-like  pigment  concentration;  b( 550) 
=  0. 1 2(pigment)0  H.  Dashed  line  Mu  represents  Morel's  upper  limit  for  total  scattering  coefficient  predicted  from 
chlorophyll-like  pigment  concentration;  6(550)  =  0.45(pigment)°  6J.  Horizontal  line  is  the  hydrosol  total  scattering 
coefficient  for  the  NOARL  blue  water  model. 


of  the  organic  detritus  component  to  the 
highest  chlorophyll  concentrations  in  the 
simulation.  The  result  of  the  curve  fit  and 
extrapolation  gave  hydrosol  total  scattering 
coefficient  values  that  were  slightly  higher 
than  Morel’s  empirical  limit  but  did  not 
exceed  the  data  limits  reported  by  Gordon 
and  Morel  (1983).  The  volume  scattering 
function  of  organic  detritus  was  estimated 
from  Petzold’s  (1972)  measurements  of  the 
volume  scattering  function  of  San  Diego 
Harbor  water  with  the  volume  scattering 
function  of  molecular  water  (Morel  1974) 
subtracted. 

Absorption  coefficients  of  the  simula¬ 
tions—  The  absorption  coefficient  for  mo¬ 
lecular  water  was  taken  from  the  tables  of 
Smith  and  Baker  (1981).  Quartz  was  as¬ 
sumed  to  be  nonabsorbing.  The  absorption 
coefficient  for  Platymonas  was  calculated 
from  Bricaud  and  Morel  (1986)  where  the 
specific  absorption  coefficient  is  normalized 


to  the  concentration  of  Chlorophyll-like 
pigments.  The  absorption  coefficient  for  or¬ 
ganic  detritus  was  estimated  from  a  regres¬ 
sion  relation  in  Prieur  and  Sathyendranath 
(1981)  which  is  normalized  to  the  total  scat¬ 
tering  coefficient  of  the  organic  detritus  at 
550  nm.  Yellow  substance  absorption  was 
estimated  by  the  assumption  that  its  ab¬ 
sorption  coefficient  was  20.7%  of  the  total 
hydrosol  absorption  coefficient  as  deter¬ 
mined  from  the  regression  analysis  of  Prieur 
and  Sathyendranath  (1981). 

NOARL  Mo.ite  Carlo  simulations 
The  NOARL  optical  model  is  a  Monte 
Carlo  simulation  of  the  radiative  transfer 
equation  (Eq.  1)  written  in  structured  For¬ 
tran  77  (Etter  1987).  The  fundamental  al¬ 
gorithm  is  based  on  the  work  of  Gordon 
and  Brown  (1973)  with  additions  from  Plass 
and  Kattawar  (1972)  and  Kirk  (1981a).  In 
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brief,  a  photon  enters  the  water  and  a  ran¬ 
dom  number  generator  chooses  whether  the 
photon  is  absorbed  or  scattered,  based  on 
the  absorption  coefficient  and  total  scatter¬ 
ing  coefficient  for  the  hydrosol.  The  random 
number  generator  subroutine  uses  a  modi¬ 
fied  random  number  shuffler  (Press  et  al. 
1986)  to  minimize  serial  correlations.  If  an 
absorption  event  is  chosen,  the  photon  is 
terminated  and  a  new  photon  enters  the 
water.  If  a  scattering  event  is  chosen,  the 
random  number  generator  is  reactivated  and 
it  selects  the  scattering  event  from  a  possible 
interaction  with  a  water  molecule,  quartz¬ 
like  particle,  algal  cell,  or  organic  detritus 
particle  depending  on  the  value  of  the  total 
scattering  coefficient  entered  for  each  class 
of  particle. 

After  the  interacting  particle  is  selected, 
the  random  number  generator  is  reactivated 
and  the  new  trajectory  of  the  photon  deter¬ 
mined  from  the  normalized  volume  scat¬ 
tering  function  (scattering  phase  function  of 
Gordon  and  Brown  1973)  for  the  interacting 
particle.  A  further  reactivation  of  the  ran¬ 
dom  number  generator  determines  the  op¬ 
tical  path  based  on  the  sum  of  the  absorp¬ 
tion  and  total  scattering  coefficient.  The 
photons  activate  counters  at  specified  depth 
intervals  intersected  by  the  photon  trajec¬ 
tory,  and  the  various  photon  sums  are  used 
to  calculate  zonal  radiances,  downwelling 
and  upwelling  irradiances,  and  downwelling 
and  upwelling  scalar  irradiances.  The  ac¬ 
curacy  of  the  photon  fluxes  is  checked  by 
calculating  the  absorption  coefficient  from 
the  recorded  fluxes  using  the  three-param¬ 
eter  model  solution  of  the  Gershun  equa¬ 
tion.  The  model  output  was  accepted  for  the 
depths  from  the  surface  down  to  the  depth 
where  the  calculated  absorption  coefficient 
deviated  by  no  more  than  2%  from  the  input 
absorption  coefficient.  Model  output  below 
this  depth  was  judged  to  have  too  high  a 
random  variance.  Duplicate  simulation  runs 
were  made  for  each  concentration  of  chlo¬ 
rophyll-like  pigment  chosen  and  each  run 
used  2.5  x  106  photons. 

The  NOARL  optical  model  was  imple¬ 
mented  as  follows.  A  flat  water  surface  was 
used  with  no  skylight  so  that  the  only  pho¬ 
tons  entering  were  of  solar  origin.  On  clear 
days  the  solar  beam  contributes  ~85%  of 


the  entering  irradiance  (Kirk  1981ft).  The 
solar  beam  was  set  at  an  11°  zenith  angle  in 
air;  the  wavelength  chosen  for  the  simula¬ 
tions  was  440  nm.  The  results  from  the  ir¬ 
radiances  generated  for  440  nm  should  not 
vary  widely  from  those  of  other  wavelengths 
in  the  400— 490-nm  waveband.  We  restrict 
consideration  of  our  simulation  without  in¬ 
ternal  radiation  sources  to  this  waveband 
because  the  500-600-nm  waveband  ex¬ 
hibits  Raman  emission  effects  in  the  clearest 
ocean  waters  (Marshall  and  Smith  1988; 
Stavn  and  Weidemann  1988a,ft),  while  flu¬ 
orescence  emission  strongly  affects  radiant 
flux  in  the  600-nm  +  region  (Gordon  1979; 
Morel  and  Prieur  1977).  Internal  source 
terms  can  always  be  added  to  simple  irra¬ 
diance  equations  after  the  situation  without 
internal  emission  has  been  properly  eluci¬ 
dated. 

The  purpose  of  this  simulation  of  the  ra¬ 
diative  transfer  equation  is  to  estimate  the 
shape  factors  defined  by  Aas  (1987)  which 
cannot  be  deduced  directly  from  simple  ir¬ 
radiance  measurements.  From  the  zonal 
(integrated  over  2ir)  radiances  at  prescribed 
depth  intervals  and  the  hydrosol  volume 
scattering  function  we  can  evaluate  the  in¬ 
tegral 


2ir  I  I  0(9,  </>,  6',  <p')  do> 

J 2wd 

■L(z\  8',  27r')sin  9'  d9', 

which  differs  from  the  numerator  of  Eq.  9 
in  having  the  radiance  L(z\  6,  <t>)  approxi¬ 
mated  with  the  zonal  radiance  I.(z\  6,  2ir). 
This  approximation  of  the  radiance  is  a 
common  method  of  solving  radiative  trans¬ 
fer  equations  in  a  less  computation-inten¬ 
sive  fashion,  and  it  is  done  both  for  empir¬ 
ical  work  and  in  the  theory  of  radiative 
transfer  (Zanevcld  1974;  Kirk  1983;  Kour- 
ganoff  1963).  The  zonal  radiance  distribu¬ 
tion  is  a  good  approximation  of  the  radiance 
distribution  established  at  solar  angles  near 
the  zenith  direction. 

The  zonal  radiance  distribution  is  a  less 
accurate  estimate  of  the  radiance  distribu¬ 
tion  established  at  greater  deviation  of  the 
solar  bearing  from  the  zenith.  However,  Kirk 
(1984),  using  a  Monte  Carlo  model  with 


Shape  factors 


1435 


i 


zonal  radiances,  has  successfully  simulated 
optical  relationships  of  various  hydrosols 
dependent  on  the  volume  scattering  func¬ 
tion  and  radiance  distribution  at  large  solar 
zenith  angles.  We  are  assuming  that  the  re¬ 
sults  obtained  in  this  study  at  a  large  solar 
zenith  angle  are  adequate  approximations 
of  results  that  would  be  obtained  with  a  full 
radiance  distribution  with  azimuth  angles. 
This  assumption  deserves  further  study. 
This  integral,  which  uses  zonal  radiances,  is 
the  backscattered  scalar  irradiance  from  the 
downwelling  scalar  irradiance.  A  similar 
term  for  the  upwelling  scalar  irradiance  re¬ 
sults  from  replacing  d  with  u  and  u  with  d. 
We  then  obtain  the  radiance  backscatter 
coefficient  for  downwelling  radiant  flux  from 
the  integral  and  the  downwelling  scalar  ir¬ 
radiance 


Bfz)  = 


P(0,  0,  6',  <j>')  da) 


■L(z;  6',  2ir')sin  6'  dfl'j 

-  E^z).  (11) 


Similarly  for  the  upwelling  radiant  flux 


bfz)  = 


0(6,  0,  O' ,  0')  do; 


■L(z;  9 ’,  2ir')sin  6'  dfl’j 
-£0u(z).  (12) 


The  shape  factors  follow 


since  the  bh  coefficient  is  evaluated  from  the 
given  parameters  of  the  Monte  Carlo  sim¬ 
ulation. 

To  answer  the  questions  posed  by  Prei- 
sendorfer  and  Mobley  ( 1 984)  and  quantified 
by  Aas  (1987),  we  have  simulated  the  ra¬ 
diative  transfer  for  the  blue,  blue-green,  and 


green  Case  1  waters  and  determined  the  r 
coefficients,  using  Eq.  1 1  and  1 2  on  the 
Monte  Carlo  output. 

Results 

Two  situations  were  modeled  to  generate 
the  suite  of  optical  water  types  in  the  Morel 
optical  classification  and  obtain  simulated 
irradiance  fields  that  may  represent  the  ex¬ 
treme  values  obtainable. 

Platymonas  new  bloom:  The  surface  lay¬ 
ers  of  a  parcel  of  blue  water  are  ventilated 
with  nutrient-rich  water  and  this  initiates  a 
bloom  of  Platymonas  cells.  Growth  is  as¬ 
sumed  rapid  enough  that  no  significant 
amounts  of  yellow  substance  or  organic  de¬ 
tritus  are  generated.  Traganza  (1969)  de¬ 
scribed  bloom  situations  in  the  Sargasso  Sea 
that  rapidly  convert  blue  waters  to  a  green 
water  parcel  and  thereafter  generate  detrital 
organic  matter.  The  pattern  of  formation  of 
green  Case  1  waters  often  seems  to  be  one 
of  a  thick  algal  layer  at  the  surface  overlay¬ 
ing  essentially  blue  water  (Austin  1 980).  This 
situation  is  also  possible  near  certain  coasts 
where  blue  water  can  receive  a  rapid  influx 
of  nutrients  from  upwellng,  etc.  This  Platy¬ 
monas  growth  is  postulated  to  be  rapid 
enough  to  increase  the  chlorophyll  concen¬ 
tration  from  neglible  amounts  to  20  mg  nr  \ 

Platymonas  old  bloom:  At  the  concen¬ 
tration  of  chlorophyll-like  pigments  of  20 
mg  m-3,  it  is  assumed  that  the  nutrients  are 
exhausted  and  the  bloom  declines.  During 
the  period  of  decline  the  decay  of  old  cells 
yields  dissolved  yellow  substance,  a  de¬ 
crease  in  chlorophyll-like  pigments,  and  an 
increase  in  organic  detritus  representing  de¬ 
composed  cells,  cellular  fragments,  and  bac¬ 
teria.  This  would  generate  the  general  con¬ 
ditions  for  Case  1  waters  postulated  by 
Morel.  The  scattering  coefficient  at  550  nm 
represents  a  higher  concentration  of  organic 
detritus  relative  to  live  cells  at  lower  con¬ 
centrations  of  chlorophyll-like  pigments  and 
a  relatively  lower  concentration  of  organic 
detritus  at  higher  pigment  concentrations. 

The  presumption  here  is  that  values  of 
optical  parameters  for  intermediate  cases 
would  fall  between  the  range  established  for 
the  Platymonas  old  bloom  and  new  bloom. 
The  new  bloom  is,  admittedly,  the  least  re¬ 
alistic  situation. 
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The  shape  factors  in  each  simulation  were 
calculated  both  for  individual  components 
of  the  hydrosol  and  the  total  hydrosol.  The 
shape  factor  for  the  total  hydrosol  would  be 
determined  from  radiance  data  obtained 
from  sensors  directly  measuring  the  radi¬ 
ances  in  the  various  optical  water  types.  The 
values  of  the  shape  factors  reported  here  are 
averages  from  the  surface  to  the  maximum 
depth  of  the  simulations  with  acceptable 
statistical  variation.  The  shape  factors  in 
this  depth  range  varied  by  only  a  few  per¬ 
cent.  Shape  factor  values  for  simulations  at 
large  solar  zenith  angles  are  not  reported  as 
they  were  only  a  few  percent  greater  than 
the  values  calculated  from  simulations  at  an 
IT  solar  zenith  angle. 

Molecular  water:  This  component  serves 
as  a  check  on  the  adequacy  of  the  model 
implementation  and  the  assumptions  used 
to  calculate  the  shape  factor.  Under  all  con¬ 
ditions  the  rd  for  molecular  water  varied 
from  1.015  to  1.022,  indicating  the  least 
accurate  estimates  of  rd  would  be  in  error 
by  2.2%  and  the  most  accurate  estimates  by 
1 .5%.  The  ru  coefficients  for  molecular  water 
varied  from  1 .002  to  1 .003  so  that  the  worst 
error  of  estimate  for  ru  coefficients  was  0.3%. 
The  greatest  error  of  estimate  with  a  solar 
zenith  angle  of  70°  was  2.3%  for  rd  and  the 
r„  error  was  unchanged.  We  thus  have  a 
measure  of  the  error  involved  in  the  com¬ 
mon  practice  of  using  zonal  radiances  when 
calculating  solutions  for  radiative  transfer 
problems.  The  greater  error  in  the  rd  coef¬ 
ficient  is,  of  course,  explained  by  the  fact 
that  the  radiances  are  more  nearly  uniform¬ 
ly  distributed  in  the  upwelling  hemisphere 
than  in  the  downwelling  hemisphere.  This 
inequality  of  the  radiance  distributions  in 
the  upper  and  lower  hemispheres  is  of  great 
importance  in  explaining  the  following  re¬ 
sults. 

Quartz:  The  shape  factors  varied  from  1 .6 
to  2.1  for  rd  and  3.9  to  4.3  for  ru  showing 
definite  effects  of  the  forward  scattering  lobe 
on  the  radiance  backscatter  coefficient  and 
thus  the  backscattered  flux  of  the  irradiance 
field. 

Organic  detritus:  The  shape  factors  varied 
from  1.6  to  2. 1  for  rd  and  4.9  to  5.0  for  ru. 
We  see  increasing  influence  of  the  forward 
scattering  in  the  ru  results  as  the  backscat- 


tering  lobe  represents  only  ~  1%  of  the  total 
scattering  coefficient  for  organic  detritus 
(Morel  and  Prieur  1977). 

Platymonas:  The  shape  factors  varied 
from  44  to  49  for  rd  and  76  to  90  for  r„. 
The  backscattering  lobe  is  only  about  0.06% 
of  the  total  scattering  coefficient  for  this  alga 
which  is  not  the  lowest  value  recorded  for 
algal  cells  (Bricaud  et  al.  1983;  Bricaud  and 
Morel  1986).  Forward  scattering  is  the 
dominant  factor  in  the  irradiance  backscat¬ 
ter  produced  by  this  suspended  particle. 

Hydrosol:  Total  rd  and  ru  coefficients  for 
each  optical  water  type  investigated  are  in 
Fig.  3.  The  rd  coefficient  varies  from  1.3  for 
blue  water  to  1 0  for  green  water  of  20  mg 
m~3  chlorophyll  concentration.  The  ru  coef¬ 
ficient  varies  from  1.8  for  blue  water  to  20 
for  green  water  of  20  mg  m  ~3  chlorophyll 
concentration.  Although  the  r  coefficients 
for  the  individual  components  of  the  hy- 
drosols  varied  by  a  maximum  of  about  23% 
for  all  the  optical  water  types,  the  variation 
of  these  coefficients  for  all  the  hydrosols  is 
an  order  of  magnitude.  At  the  lowest  par¬ 
ticle  concentrations  the  shape  factor  is  pri¬ 
marily  an  interaction  between  the  volume 
scattering  functions  of  the  water  molecule 
and  suspended  quartzlike  particles.  Adding 
suspended  material  such  as  algae  or  organic 
detritus,  which  is  more  strongly  forward 
scattering,  increases  the  magnitude  of  the 
shape  factor  for  the  hydrosol  (Figs.  2  and  3, 
Table  1 ).  The  shape  factors  for  the  hydrosols 
show  a  tendency  to  increase  to  a  limit  im¬ 
posed  by  the  maximum  shape  factor  for  the 
dominant  added  component.  The  overall 
increase  in  the  shape  factor  of  the  hydrosols 
represents  the  increase  in  particulate  matter 
and  the  forward  scattering  of  this  material. 
This  is  correlated  with  an  increase  in  the 
hydrosol  backscatter  coefficient— the  pa¬ 
rameter  to  which  the  shape  factor  is  nor¬ 
malized. 

Discussion 

The  question  is  under  what  conditions 
can  simple  irradiance  models  allow  the  in¬ 
version  of  irradiance  data  to  yield  the  in¬ 
herent  optical  properties  of  natural  hydro¬ 
sols?  The  simplest  irradiance  model— the 
three-parameter  model  — allows  inversion 
of  the  net  irradiant  flux  to  yield  the  hydrosol 
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Fig.  3.  Downwelling  and  upwelling  shape  factors  plotted  against  pigment  concentration  of  various  ocean 
hydrosols.  Solid  lines  represent  shape  factors  for  the  Platymonas  new  bloom  and  dashed  lines  represent  those 
for  the  Platymonas  old  bloom. 


absorption  coefficient  in  any  hydrosol.  In¬ 
version  for  the  backscatter  coefficient,  how¬ 
ever,  is  problematical.  After  preliminary 
calculations  indicated  that  the  shape  factors 
were  too  high  in  turbid  Case  2  waters  to 
allow  inversion  of  the  two-flow  model,  Aas 
(1987)  noted  that  in  clear  ocean  waters  (blue 
and  blue-green  waters)  the  backscatter  from 
molecular  water  is  dominant,  which  implies 
that  the  shape  factors  might  approach  1.0 
and  the  two-flow  model  could  be  inverted. 
We  have  shown  here  that  the  smallest  val¬ 
ues  for  the  shape  factors,  obtained  from  blue 
water  (molecular  water  plus  quartz),  are  1.3 
for  rd  and  1 .8  for  ru.  The  minimum  error 
for  an  attempt  to  invert  the  two-flow  equa¬ 
tions  in  the  clearest  possible  ocean  water 
would  be  55%.  The  inclusion  of  any  other 
scattering  material  into  the  ocean  will  in¬ 
crease  the  error  dramatically  to  an  order  of 
magnitude  in  Case  1  waters  of  20  mg  m~5 
chlorophyll  concentration.  On  the  basis  of 
this  evidence  we  conclude  that  the  mean 
backscatter  coefficient  (5)  is  an  apparent  op¬ 
tical  property  that  cannot  be  equated  with 
the  inherent  hydrosol  backscatter  coeffi¬ 


cient.  Therefore  inversion  of  the  two-flow 
model  produces  a  mean  backscatter  coeffi  - 
cient  that  is  a  backscatter  coefficient  for  the 
irradiances  and  is  net  a  backscatter  coeffi¬ 
cient  for  the  hydrosol  in  which  the  irradi¬ 
ances  are  measured. 

A  more  serious  problem  arises,  however, 
when  we  consider  that  the  ratio  of  the  shape 
factors  (rjrd)  for  all  simulations  of  blue- 
green  and  green  waters  ranges  from  1.39  to 
2. 1 7  as  shown  in  Table  2.  Thus  the  radiance 
distribution  :n  natural  hydrosols,  already 
noted  as  being  rather  different  in  the  down- 
welling  and  upwelling  hemispheres,  renders 
bd  unequal  to  bu  with  the  hypothesis  of  a 
single  mean  backscattering  coefficient  un¬ 
supported.  This  extends  the  unknown  pa¬ 
rameters  of  Fq.  3  and  4  in  natural  hydrosols 
to  three  in  each  equation.  There  is,  of  course, 
no  unique  solution  for  such  a  system.  We 
tried  to  correct  for  this  by  adding  the  three- 
parameter  model  solution  to  the  two-flow 
model  so  that  we  would  have  a  system  of 
three  equations  in  three  unknowns  (Eq.  1 3- 
1 5).  A  simple  application  of  Cramer’s  rule 
for  the  solution  of  nonhomogeneous  linear 


Table  1.  Optical  coefficients  (m  ')  used  in  simulations. 
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Table  2.  Shape  factor  ratios. 


Pigment  concn 
(mg  m  ') 

Plalymonas 
new  bloom 
(fA) 

Plalymonas 
old  bloom 
('■Jr 4) 

0.0 

1.36 

1.36 

0.1 

1.39 

1.79 

0.5 

1.50 

2.03 

1.0 

1.57 

2.10 

5.0 

1.80 

2.17 

10.0 

1.93 

2.03 

15.0 

2.01 

1.98 

20.0 

1.97 

2.07 

equations  for  the  unknown  coefficients  a, 
Bd,  and  Bu  revealed  that  such  a  system  is 
inconsistent  (i.e.  has  no  solution  at  all): 


and 


dE(z) 

dz 


a(z) 

M~) 


E(z), 


dEJz) 

dz 


a(z) 

BA=) 


EA=) 


hJH 


EJ(z) 


Bu(z) 

BA?) 


Eu(z) 


AEu(z) 

dz 


g(r) 

BAz) 


EAz) 


bAz) 

BAz) 


EAz) 


+ 


ba= ) 

BAz) 


EAz). 


(13) 


(14) 


(15) 


In  this  connection  it  is  worthwhile  to  re¬ 
view  some  results  of  atmospheric  work  uti¬ 
lizing  Schuster-type  equations  for  an  at¬ 
mospheric  two-flow  model.  Kondratyev 
(1969)  reported  that  the  atmospheric  two- 
flow  model  with  a  uniform  volume  scatter¬ 
ing  function,  which  would  give  rd  =  ru  = 
1.0,  provided  an  adequate  solution.  When 
he  attempted  a  solution  in  the  atmosphere 
for  an  asymmetric  volume  scattering  func- 
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rough  estimates  of  the  backscatter  coeffi¬ 
cient  and  not  an  accurate  solution. 

It  is  possible  to  make  algebraic  substitu¬ 
tions,  however,  in  the  above  system  to  con¬ 
vert  it  to  a  form  that  does  have  a  restricted 
solution  dependent  on  the  ratio  of  the  shape 
factors— the  same  factors  which  originally 
rendered  the  equations  inconsistent— 

5-w-  **>  w, ■■ 

and  the  following  equation  system  does  have 
a  solution  if  rjrd  is  known: 


dE(z) 
d  z 


a(r) 

£(z) 


E(z), 


(16) 


d EJ,z)  a(z) 


~  bAZ\uz) 


_  ru(z)  fl(z)\ 
rj^z)  fiu(z)) 


dEu(z)  a(z) 


Bu(Z{mz)  rt(z)  UW(z)}Eu{:) 


where  R(z)  =  E„(z)/EA=) 


Equations  16-18  constitute  what  we  term 
the  net  flow-two-flow  model,  which  can  be 
inverted  to  provide  the  inherent  absorption 
coefficient  for  a  hydrosol  and  the  radiance 
backscatter  coefficients  for  the  light  field 
which  are  one  step  removed  from  the  hy¬ 
drosol  backscatter  coefficient  when  ru  lrd  is 
known.  The  model  represents  three  separate 
integrations  of  the  radiative  transfer  equa¬ 
tion  (Eq.  1).  The  iterative  techniques  de¬ 
veloped  by  Preisendorfer  and  Mobley  ( 1 984) 
appear  to  be  the  preferred  mode  for  devel¬ 
oping  the  inversion  algorithm  for  this  equa¬ 
tion  system. 

Where  do  we  stand,  then,  on  the  question 
of  inverting  the  irradiances  of  natural  light 
fields  for  the  backscatter  coefficient  of  a  hy¬ 
drosol?  The  problem  stems  from  the  asym¬ 
metry  of  the  volume  scattering  function  in 
natural  hydrosols  and  the  differences  be¬ 


tween  the  radiance  distributions  in  the 
downwelling  and  upwelling  hemispheres. 
These  aspects  of  natural  hydrosols  require 
the  use  of  shape  factors  when  inversions  are 
attempted,  and  we  have  shown  that  these 
parameters  vary  as  a  function  of  the  con¬ 
centration  of  suspensoids.  Coupling  the 
three-parameter  model  with  the  two-flow 
model  (net  ftow-two-ftow  model)  should  al¬ 
low  the  inversion  for  the  radiance  back¬ 
scatter  coefficient  in  most  green  Case  1 
waters  where  the  shape  factor  ratio  ( rjrd ) 
is  close  to  2.0  (Table  2).  There  is  a  regular 
increase  in  the  magnitude  of  the  radiance 
backscatter  coefficient  caused  by  the  in¬ 
crease  of  the  shape  factors  with  the  increase 
of  the  hydrosol  backscatter  coefficient  (Figs. 
2,  3).  This  fact  will  allow  a  functional  re¬ 
lationship  to  be  established  between  the  ra¬ 
diance  backscatter  coefficient  and  the  hy¬ 
drosol  backscatter  coefficient  from  the 
results  of  NOARL  Monte  Carlo  simula¬ 
tions.  In  this  application,  an  improvement 
of  up  to  an  order  of  magnitude  in  the  esti¬ 
mate  of  the  hydrosol  backscatter  coefficient 
is  to  be  expected  over  earlier  estimates. 

The  situation  is  not  as  simple  in  the  blue 
and  blue-green  Case  1  waters  where  the 
shape  factor  ratio  varies  downward  some¬ 
what  from  2.0,  depending  on  the  relative 
concentrations  of  algal  cells  and  organic  de¬ 
tritus.  Another  approach  to  the  inversion 
problem  has  been  proposed  by  Zaneveld 
(1982);  he  rejects  the  two-flow  approach  in 
favor  of  integrations  of  the  radiative  trans¬ 
fer  equation  that  give  a  statement  of  the 
optical  factors  affecting  the  transmittance 
and  absorptance  of  the  upwelling  nadir  ra¬ 
diance  ratioed  to  the  downwelling  scalar  ir- 
radiance.  Shape  factors  are  also  involved 
with  these  integrations  but  they  are  defined 
differently  from  Aas’  shape  factors  for  the 
downwelling  and  upwelling  irradiance 
streams.  The  solution  of  Zaneveld’s  model 
requires  the  same  up-  and  downwelling  ir¬ 
radiances  required  in  the  net  flow-two-flow 
model  plus  measurements  of  the  upwelling 
radiance  received  from  the  nadir  viewing 
angle.  The  absorption  coefficient  can  be  ob¬ 
tained  from  the  three-parameter  model  so¬ 
lution  of  the  Gershun  equation  as  is  done 
in  the  net  flow-two-flow  approach.  The  as¬ 
sumptions  required  to  solve  the  Zaneveld 
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model  are  that  the  shape  factor  for  the  back- 
scatter  of  downwelling  scalar  irradiance  to 
the  nadir  radiance  (fb)  and  the  radiance  shape 
factor  for  the  upwelling  radiances  (fL )  both 
be  equal  to  1 .0.  If  these  assumptions  are  met 
the  model  can  be  inverted  for  the  back- 
scatter  coefficient  of  the  hydrosol.  Thus,  in¬ 
version  for  the  hydrosol  backscatter  coeffi¬ 
cient  can  be  accomplished  from  the 
integrations  of  the  radiative  transfer  equa¬ 
tion  that  yield  the  net  flow-two-flow  model 
or  the  integrations  that  yield  the  Zaneveld 
model. 

The  assumptions  required  to  solve  the  net 
flow-two-flow  model  in  the  green  Case  1 
waters  and  the  assumptions  required  to  solve 
the  Zaneveld  model  both  imply  errors  in 
the  range  of  ~  30%  with  Zaneveld  claiming 
this  error  range  for  all  oceanic  waters.  How¬ 
ever,  this  error  range  was  estimated  with  a 
limited  number  of  cases  in  which  variations 
in  suspensoid  concentration  were  not  ac¬ 
counted  for  and  the  volume  scattering  func¬ 
tion  for  algal  cells  was  not  tested.  It  is  im¬ 
portant  because  the  strongest  effect  on  the 
shape  factors  determined  in  the  NOARL 
simulations  for  the  net  flow-two-flow  mod¬ 
el  was  due  to  the  volume  scattering  function 
for  algal  cells.  An  inversion  with  the  net 
flow-two-flow  model  requires  fewer  light 
field  measurements  while  an  inversion  with 
the  Zaneveld  model  requires  more  light  field 
measurements  but  may  provide  greater  ac¬ 
curacy.  NOARL  simulations  currently  in 
progress  which  have  a  greater  range  of  op¬ 
tical  parameters,  varying  solar  and  sky  con¬ 
ditions,  and  coverage  of  the  400-490-nm 
waveband  will  provide  the  data  needed  to 
investigate  variance  in  all  the  shape  factors 
proposed  so  far  and  determine  the  best  in¬ 
version  models  for  the  various  optical  water 
types. 

Conclusions 

The  two-flow  model  cannot  be  inverted 
to  obtain  the  true  backscatter  coefficient  bb 
of  natural  hydrosols.  The  net  flow-two-flow 
model  at  present  can  be  inverted  in  all  hy¬ 
drosols  to  obtain  the  absorption  coefficient 
a.  When  the  ratio  of  the  shape  factors  rjrd 
is  known,  it  has  the  potential  to  be  inverted 
for  the  radiance  backscatter  coefficients  Bd 
and  Br  The  NOARL  simulation  of  the  ra¬ 


diative  transfer  equation  provides  estimates 
of  the  shape  factors  rd  and  ru.  The  shape 
factors  ./J,  and  fL  can  also  be  estimated.  The 
estimates  of  rd  and  r„  will  establish  the  basis 
for  a  functional  relationship  between  the  hy¬ 
drosol  backscatter  coefficient  bb  and  the  ra¬ 
diance  backscatter  coefficients. 
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